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The quantum corrections to the conductivity and the thermopower in monolayer graphene are 
studied. We use the recursive Green's function method to calculate numerically the conductivity 
and the thermopower of graphene. We then analyze these weak localization corrections by fitting 
with the analytical theory as function of the impurity parameters and the gate potential. As a result 
of the quantum corrections to the thermopower, we find large magnetothermopower which is shown 
to provide a very sensitive measure of the size and strength of the impurities. We compare these 
analytical and numerical results with existing experimental measurements of magnetoconductance 
of single layer graphene and find that the average size and strength of the impurities in these 
samples can thereby be determined accurately. We suggest favorable parameter ranges for future 
measurements of the magnetothermopower. 

PACS numbers: 65.80.Ck, 72.20.Pa, 73.22. Pr 
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I. INTRODUCTION 



Since its synthesis, graphene^ has attracted a lot of 
attention, both for its novel electronic properties and its 
possible applications. One of the most remarkable as- 
pects is that the weak localization effect, which is a typ- 
ical low-temperature phenomenon due to quantum cor- 
rections to the conductivity 5a, can in graphene be ob- 
served up to a temperature range of 200 Kelvin. ^^"^ It has 
been theoretically predicted that the sign of these cor- 
rections strongly depends on the kind of impurities in 
the graphene sample,^ so that both positive and nega- 
tive magnetoconductivity. the so-called weak antilocal- 
ization effect, can be observed. This effect is well known 
to occur due to spin-orbit interaction, even though the 
latter is typically very weak in graphene. Noting that the 
graphene lattice is composed of two sublattices, one can 
formulate the sublattice degree of freedom as an isospin 
which is strongly coupled to the momentum. Thereby, 
any elastic scattering that breaks the graphene sublat- 
tice symmetry is expected to result in weak antilocaliza- 
tion. This is made more complicated, however, by the 
fact that there is another discrete degree of freedom in 
graphene, the two degenerate Dirac cones, which can cor- 
respondingly be formulated by introducing a pseudospin 
index. Accordingly, any scattering which mixes these two 
valleys can result in yet another change of the sign of 
the quantum correction to the conductivity, and thereby 
to the restoration of weak localization. More recently, 
thermopower S in graphene has been measured by sev- 
eral groups.^"* The classical thermopower of graphene 
has been calculated with a range of different methods. 
The experiments show good agreement with the Mott 
formula, which corresponds to the leading term of the 



Sommerfeld expansion at low temperatures and large 
gate voltages.^ Close to the Dirac point, thermopower 
shows an unusual behavior, being linear in gate voltage, 
and changing sign. Thus, while usually thermopower 
is expected to increase as the Fermi energy, the elec- 
tron density and thereby the conductivity are lowered, in 
graphene the thermopower becomes smaller as the Dirac 
point is approached and the electron density and the con- 
ductivity are reduced. 

The main aims of this paper are, firstly, to find out 
if there are sizable quantum corrections to the ther- 
mopower SS in graphene and whether they are sensitive 
to weak magnetic fields, secondly, to present and analyze 
the combined analytical and numerical theory in order to 
provide the basis for a more quantitative analysis of ex- 
perimental results on magnetoconductance of graphene. 
This will allow one to characterize samples, in particular 
the density, strength, and size of carrier scatterers, more 
accurately. 

We restrict our attention to the diffusion term of the 
thermopower and do not consider the phonon drag con- 
tribution, which becomes relevant only at high tempera- 
tures. In normal metals, the quantum corrections to the 
diffusion thermopower are known to be dominated by the 
weak localization corrections to the conductivity, yield- 
ing SS/S sa —Sa/a}'^^^'^ Thus, when the conductivity 
is reduced by the weak localization correction, the ther- 
mopower becomes enhanced. This is expected since the 
thermopower is known to increase from a metal towards 
an insulator. Because of the particular electronic proper- 
ties of graphene, where the sign of the quantum correc- 
tion can change upon varying the gate voltage, it could be 
expected that these quantum corrections to thermopower 
are particularly large in graphene. We therefore deter- 



mine the strength and the sign of the resuhing magne- 
tothermopower and analyze in detail how it changes with 
impurity parameters such as range, density, and strength 
and with the gate voltage. 

This paper is organized as follows. We start in Sec. 
II with a short introduction to the electronic properties 
of graphene. We introduce the model for the description 
of the impurity potential in Sec. Ill and list the result- 
ing momentum scattering rates and their dependence on 
the impurity matrix elements. In Sec. IV A we give a 
brief review of transport theory, in particular the the- 
ory of weak localization corrections to the conductivity 
and its application to graphene. In Sec. IV B we review 
the theory of thermopower and address leading quantum 
corrections. 

In Sec. V we present a numerical method to calculate 
the conductance which is based on the recursive Green's 
function technique (Sec.VB). We introduce the Hamilto- 
nian that models impurities in graphene for the numerical 
calculation in Sec. VC and relate it in Sec. VI to the 
impurity model introduced in Sec. III. This connection 
makes it possible to display the related transport scat- 
tering rates from Sec. Ill and Sec. IV as functions of 
the sample parameter of the numerical method. Sec. V, 
which is done in Sec. VI B. 

In Sec. VII we present the numerical results for the 
conductance, Sec. VII A, and thermopower. Sec. VII B. In 
Sec. VIII we analyze the numerical results by fitting 
them to the analytical results in Sec. VIII A. In Sec. 
VIII B we attempt an ab initio analysis. To this end, we 
use the relation of the scattering rates to the impurity 
parameters of the numerical calculations and insert it in 
the analytical formulas given in Sec. IV. In Sec. VIII C 
we present the analytical results for the quantum correc- 
tions to the conductance at zero magnetic field. In Sec. 
VIII D we present the results for the quantum correc- 
tions to the thermopower at zero magnetic field and the 
resulting magnetic field dependence of the thermopower. 
In Sec. IX we compare the numerical and the analyt- 
ical results with experimental results on weak localiza- 
tion corrections to the conductance. Finally, we draw 
the conclusions and summarize our results in Sec. X. 
In appendix A the Hamiltonian is given in matrix nota- 
tion. In appendix B we present analytical results for the 
magnetoconductance when the warping rate is neglected, 
1/t,„ = 0. 



II. ELECTRONIC STRUCTURE OF 
GRAPHENE 

Let us start with a brief review of the electronic 
properties of graphene, introducing the notation used in 
this article. We follow the convention of McCann and 
coworkers.^ For an overview and comparison of notations, 
see also 13, 14 and 15. 

Graphene is a two-dimensional layer of carbon atoms 
arranged in a honeycomb lattice. These carbon atoms are 



connected by strong cr-bonds with their three neighboring 
atoms. The corresponding energy bands are filled valence 
bonds, lying deep below the Fermi energy. The 7r-bond 
leads to the formation of a tt— band which is exactly half 
filled in ungated and undoped graphene. Thus, to study 
transport properties we can restrict the model to this tt- 
band. The basis of the Bravais lattice consists of two 
atoms, each one of them forming one sublattice, named 
A- and B-sublattice. 

The electronic band structure shows two degenerate 
half filled cones, which in momentum space are located 
on the sites of the hexagonal reciprocal lattice. Close to 
zero energy, the Dirac points, the dispersion is in good 
approximation linear. Each Dirac cone is part of a sub- 
lattice, the K- and K'-valleys. 

Accordingly, the Hamiltonian can be written as a 4 x 4- 
matrix which, close to the Dirac points in linear approx- 
imation, becomes 



Hi = vfCs ® ct • fc. 



(1) 



where C3 is the diagonal Pauli matrix in K-K' space, k^ 
and ky are components of the momentum k in the plane, 
and 



(2) 



with (71.2 the non-diagonal Pauli matrices in the sublat- 
tice (A,B)-space. The next correction term to the dis- 
persion is quadratic and given by 

H2 = -M h (kl - kl) - 2(72 (k^ky)] Co, (3) 

where ^ is the parameter of warping, which is given 

byl6,17 



M 



3iag 



(4) 



In this notation the 4-component Bloch-states are given 
by 



^■' = ($ 



AK, ^BK, ^BK 



,'^A 



K' 



(5) 



The Hamiltonian can be more compactly written in the 
isospin and pseudospin notation, which is particularly 
convenient for the purpose of understanding the sign of 
weak localization corrections to the conductivity. Follow- 
ing the definition of McCann and coauthors,^ we set 



and 



where Ui, i = 1,2,3 are the Pauli matrices and cto the 
identity matrix in the sublattice (A,B)-space, and Q, i = 



^0 = Co ® CTo, 


^1= Cs ® CTi, 


(6) 


^2 = Ca ® 0-2, 


^3= Coders, 


(7) 


Ao = Co ® CTo, 


Ai==Cl®cr3, 


(8) 


A2 = C2 ® 0-3, 


A3= Cs ® o-Q, 


(9) 



1,2,3 are the Pauli matrices and Co the identity matrix 
in K-K' space. The vector of the 4x4 matrices S^ = 
(Si,S2,S3) is the zsos;jm vector referring to the suWattice 
(A.B) degrees of freedom, and A-^ = (Ai, A2, A3) is the 
pseudospin vector referring to the valley (K, K') degrees 
of freedom. 

In this notation the Hamiltonian Hi + H2 can be written 
as 




Jirac cone 



Warping correction to the Dirac cone 



(10) 



Note that Hi is independent of the pseudospin A , while 
the warping term H2 breaks the pseudospin symmetry. 



where we introduced _L = 1,2 to denote the transverse 
components, noting that r^^ — t^ 



T23 ==: ^j_j > since V{j 



1/2 

"^2 J 



T^^ and T^/ 
Vf , and V?, = V^2 



V' 



The intravalley scattering rate is accordingly given by 

(15) 



r^' ■■= 2rl] 



'33 • 



The trigonal warping term in the kinetic energy of 
graphene results in an asymmetry of the dispersion at 
each valley with respect to momentum inversion,^® while 
the total Hamiltonian including both valleys preserves 
that symmetry. Therefore, the pseudospin is expected 
to precess in the presence of the warping term. In the 
presence of elastic scattering the pseudospin channels are 
expected to relax accordingly, as in motional narrowing, 
in proportion to the scattering time r. 



III. IMPURITIES IN GRAPHENE 

Disorder can have many different origins in graphene 
such as adatoms and vacancies, while substitutional de- 
fects are rather unlikely due to the strength of the a- 
bonding. For the purpose of the transport calculations, 
it is rather convenient to classify the disorder according 
to whether it breaks the isospin symmetry (A, B sublat- 
tices) or the the pseudospin symmetry (K, K' valleys). 
Thus, the Hamiltonian of nonmagnetic disorder has in 
the representation of the pseudospin A and the isospin 
S the general form,'* 



Hi, 



/K),o(r) + 






s.Aj-y,,, {T 



(11) 



where / is the identity matrix. Vg^o (^) is the part of the 
disorder potential which leaves both the isospin and the 
pseudospin symmetry invariant, while the other terms 
are breaking either of these symmetries with the ampli- 
tudes Vij (7). We note that this disorder Hamiltonian 
is invariant under time reversal, A — >■ — A and S — >■ — S 
and thus indeed describes nonmagnetic disorder. See Ap- 
pendix A for the explicit matrix representation of the 
Hamiltonian. 

The corresponding scattering rates are given by 



(^^^.j) /^' 



(12) 



where v is the density of states per spin and valley. We 
can thus define the total scattering rate as 



i,j=l,2,3 



(13) 



The intervalley scattering rate is obtained by summing 
over all matrix elements which couple to the transverse 
pseudo spin components Ai and A2, yielding 






2r^^ 



(14) 



Tji = 2r 






(16) 



IV. 



TRANSPORT THEORY 



A. Weak Localization Corrections to the 
Conductance in Graphene 

Quantum corrections to the conductance, the so-called 
weak localization corrections, originate from the quan- 
tum interference of electrons propagating through the 
sample. On time scales t exceeding the elastic mean free 
time T, electrons move diffusively and their return prob- 
ability can be enhanced due to the constructive interfer- 
ence of the amplitudes for propagation on closed diffusion 
paths, as long as that time t does not exceed the phase 
coherence time t^{T). In the presence of spin-orbit inter- 
action, the spin precesses as the electron moves on closed 
paths, and the phase of the amplitudes for clockwise and 
anticlockwise propagation no longer match, unless the to- 
tal spin of these two propagations adds up to zero. This 
spin singlet channel leads to destructive interference, and 
thereby an enhancement of the conductance, the weak 
antilocalization correction. 

Since the spin-orbit interaction in graphene is weak, in 
the absence of magnetic impurities, there are two inde- 
pendent spin channels in graphene and the quantum cor- 
rections are doubled due to this spin degeneracy. How- 
ever, since the electron momentum jf is directly coupled 
to the isospin S (due to the A-B sublattice degree of 
freedom) in the electron Hamiltonian of Eq. (10), any 
momentum scattering will result in the breaking of the 
isospin symmetry. Therefore, there is a finite contribu- 
tion from the interference of two closed electron paths 
whose total isospin is zero, leading to weak antilocaliza- 
tion in the isospin singlet channel. 

Let us now consider the influence of the pseudospin 
degree of freedom of the two valleys in graphene on the 



weak localization corrections. Formulating the conduc- 
tance corrections in the representation of pseudospin sin- 
glet and triplet modes, one has, in two dimensions,^ 



Sa^ 



2e^De f (fQ 



nh 



(2^) 



(-C° + Ci + C2 + Co^) , (17) 



where D^ is the diffusion constant, which in graphene is 
related to the elastic scattering time as D^ = VpT, and 
the momentum integral has an upper cutoff l//e, where 
le = VpT is the elastic mean free path. The superscript 
corresponds to the pseudospin and the subscript to the 
isospin. 

Without magnetic field {B = 0) the Cooperon modes 
Cq are given by 



Thus, the pseudospin triplet Cooperon modes are at- 
tenuated. Performing the integral over momentum one 
thus obtains logarithmic weak localization corrections at 
B = 0, 



where 



6a 



T,f, T0 T^ 
Ttr ' Ti' T^ 



—0 ( ]± Zl I± 

Trh I Ttr ' Ti' T^, 



2 In 



r0 



In 1 



Ttr 

-2 In II 



J'P 



1± + I± 



(23) 



(24) 



c'oiQ) 



1 



[d,q^ 



n 



(18) 



Here, tT (T) is the dephasing rate caused by electron- 
electron and electron-phonon interactions, which pro- 
vides the low-energy cutoff of the diffusion pole of Eq. 
(18). The elastic scattering from impurities can break the 
pseudospin symmetry and results in the following pseu- 
dospin relaxation rates, 



where 



0, 

2t; 



^ — ^ — ~i 






(19) 
(20) 
(21) 

(22) 



and Ttr ~ 2r. The dependence of the various scattering 
rates [Eqs.(14). (15) and (16)] on the impurity potential 
amplitudes Vij is given by Eq. (12). One can see that 
both the sign and the amplitude of the weak localization 
corrections depend strongly on the impurity type. 

As we will study in detail below, impurities with large 
correlation lengths mix valleys weakly, l/r^ w 0, and 
therefore only attenuate two of the pseudospin triplet 
modes, in a similar way as the relaxation rate due to 
the warping term of Eq. (16). This effect leads to the 
vanishing of the weak localization effect and a flat mag- 
netoconductance. 

Upon applying an external magnetic field, the weak 
localization corrections are suppressed. Solving the 
Cooperon equation in the presence of a magnetic field by 
summing over the Landau levels, one finds that Aa{B) = 
<j{B) — a{B = 0) is given by 



Aa(i?) = — 
tt/i 



Pseudospin 
singlet 



Pseudospin triplet 




where the function F[z) in Eq. (25) is given by 

Ffuu(^) = -V' (2 + ^) +'^' (^ + l) +ln[2rB/rt,], 

(26) 
where ^'(•2) is the Digamma function and the magnetic 
rate is 



AeD.B 



For weak magnetic fields, one can use the simplified form'^ 



F(z) = ln(z) + 7/. (1/2 + z^i) 



Equation (25) can be expanded for small magnetic 
fields, yielding 



Afj(S) 
where 



.«A'^^^)^fe-^). (»> 



(27) /3(^,^ 



T* n 



1 + 2^ 



\ -i- lA -1- lA 

Ti T^ 



(30) 



(28) In Fig. 1 we plot the curve obtained from the solution 
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Figure 1. t^/t^ - r^/ri diagram illustrating the transition 
from positive magneto conductance (PMC) to negative mag- 
neto conductance (NMC) as obtained from Eq. (30). The 
transition is indicated by the dashed black line, correspond- 
ing to Aa{B) = 0. This behavior is in good agreement with 
experiments."^ 



of A(j{B) = as function of the parameters t^/t^. and 
Ttj,/Ti. This curve separates the region of positive and 
negative magnetoconductance. We note that this separa- 
tion line does not coincide with the curve S(j{B = 0) = 0, 
which corresponds to the vanishing of the weak local- 
ization correction and separates the parameter space re- 
gions of weak localization and weak antilocalization (cor- 
responding to negative and positive quantum corrections 
to the conductivity at i? = 0, respectively). As seen in 
Eq. (23) the latter separation line depends on an addi- 
tional parameter, the total scattering rate l/r. 



klT 



d\n{a{E)) 
dE 



(31) 



E=Ef 



This formula is valid for low temperatures, T <C Ep and 
large chemical potential. Here, e = — |e| is the negative 
electron charge. Thus, when the carriers have negative 
charge, thermopower is expected to be negative, while 
for holes it becomes positive. We note that the elemental 
unit of thermopower is given by a ratio of natural con- 
stants 5*0 = /cB/|e| ~ ^Q^V/ K . The thermopower can be 
also written in terms of the back gate voltage Vbg as 



S 



fc|ri da dVBG 



a dVBG dE 



(32) 



E — Ep 



In monolayer graphene, the thermopower is well de- 
scribed by Eq. (32) in a wide range of temperatures 
and gate voltages not too close to the Dirac point. ^ We 
note that this formula is valid not only classically but 
also includes quantum corrections through the conduc- 
tivity. Thus, expanding in the quantum correction to 
the conductivity ba. we can write the leading quantum 
corrections to the thermopower, (JS*, as 



bS_ 



da 
a 



dda{E) 



dE 



I 



E=Ef 



da{E) 
dE 



(33) 



E=Ef 



This relation is also supported by direct diagrammatic 
calculations of quantum corrections. ^°"^^ In standard 
metals, the last term on the r.h.s of Eq. (33) is small 
and the thermopower is dominated by the quantum cor- 
rection to the conductivity. In this work we revisit this 
relation to find out whether this also holds for graphene, 
or whether the second term in Eq. (33) is sizable. 



THE RECURSIVE GREEN'S FUNCTION 
METHOD 



B. Thermopower 

Applying a thermal gradient VT to a metallic sample 
induces not only a thermal current but also an electri- 
cal current density j = a{E + V/i/e) — rfs/T, where 
E is an applied electric field and (x is the chemical po- 
tential. Here, ?? denotes the thermoelectric coefRcient. 
Under open circuit conditions, the thermal gradient re- 
sults in a finite voltage U proportional to the temperature 
difference AT, where the proportionality constant is the 
thermopower: S ~ f^/c. At low temperatures, the ther- 
mopower is dominated by the diffusion of the electrons 
in the sample and the phonon drag contribution becomes 
small. Expanding in the ratio of temperature to Fermi 
energy and keeping the leading term, one arrives at the 
Mott formula, which relates the thermopower S to the 
derivative of the conductivity with respect to the Fermi 



In this section we introduce the numerical method 
used to calculate the electrical conductivity and the 
thermopower. We employ the recursive Green's func- 
tion method, ^^'^*', which has been previously applied to 
graphene in Refs. [21-24]. Below, we begin by reviewing 
the essential elements of this method. 

The graphene sample is assumed to be connected to 
two semi-infinite leads which are modeled by a square 
lattice. When contacting the leads at the zigzag edges of 
the graphene sample, there is no wave function mismatch 
between the propagating modes in the square lattice leads 
and the graphene sample, provided that a proper energy 
shift is used.^'^ The graphene sample is sliced transversely 
into N equal cells, with each cell containing M sites. We 
study the transport as function of N by changing the 
length L of the sample, and as function of M by changing 
the sample width W . When the free edges are of the 
armchair type, the sample dimensions are related to to N 



and M by 



V3 (N_ 
2 V 2 



1 



6 






M — 1. where 



ao = 2A6A is the lattice constant set as the distance 
between two atoms of the same sublattice. See Fig. 2 . 



A. The conductivity 

In two dimensions, the electrical conductivity a is re- 
lated to the conductance G by the standard expression 
a — LG/W and thus can be expressed in terms of re- 
tarded Green's function through the Caroli formula 



£2e2 
W~h 



-Ti-s (riG^r^G^ 



(34) 



where G^'^ denotes the advanced/retarded Green's 
function connecting two opposite contact regions (right 
and left). Here, Tr^ denotes the trace over transverse 
sites at the lead-sample interface. The matrix elements 
of the level width matrix Tp can be expressed in term of 
the transverse wave functions Xi'(*) of the lead propagat- 
ing modes. 



rp(«,«') = y'x^(0— ^Xi'(»')- 



(35) 



Here, the sum runs over the propagation modes v and v^, 
is the longitudinal propagation velocity. Next we describe 
the method used to obtain the Green's functions of the 
sample. 



B. The Green's function 

In order to obtain the Green's function amplitude 
G^ (j,j) between sites i,j at the contacts p,q, we start 
with the surface Green' function of one of the contacts, 
which is presumed to be known. Then we add one slice 
of transverse sites i = 1, .., M, and calculate the Green's 
function from the contact up to that slice. We repeat the 
procedure, adding slice by slice, until the end of the sam- 
ple is reached. This procedure is done from left to right, 
and then in the opposite direction in order to calculate 
the full Green's function of the system. The slicing is 
displayed in Fig. 2. 

Each line in Fig. 2 represents a hopping amplitude 
from one sublattice to another. This is described by 
the hopping matrix U. This matrix, together with the 
Green's function of every new single slice f/„ and the 
Green's function that includes the contacts and the sam- 
ple up to slice n, defining G'*^', are inserted in the Dyson 
equation G = G^°^ + G'^^^VG to obtain the new G. One 
can easily understand that the speed of the procedure 
strongly depends on the slicing scheme adopted. In gen- 
eral, more slices mean more steps to calculate, but each 
slice consists of fewer atom sites and so its calculation is 
done faster since an inversion at each slice is required. 
The complexity of the calculation scales as 0{NM^). 



The Green's function of the contacts has been derived 
previously and is given by'^^ 



{2\t^\f 




with 



Pu = -Vgato + 2|f^|cos 



Af +1 



iO^ 



(36) 



(37) 



where v = 1, ...,M represents the modes of the transverse 
wave function and i* is the hopping rate in i-direction. 
We can transform G''°""(j^) from the channel representa- 
tion into G''''""(j, j') in the site representation by 



M 



G^""(j,/) = Ex^(j')G"™(^)x-(j")' 



u=l 



with the transverse wave functions Xv{j) given by 



Xuij) 



-sm 



A/ + 1 V^-^ + 1 



(38) 



(39) 



C. Inclusion of impurities 

We now look at the model disorder used in the nu- 
merical calculations. The sites of the graphene sample 
are occupied by uniformly randomly distributed Gaus- 
sian scatterers with a random potential Vn G [—51^,51^]. 
We can write for the overall potential resulting from all 
-^imp scatterers at the points i?„ as 



Ni, 



vi^)^ E^»(^ 



J2Vne 



k7-«^r 



(40) 
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Figure 2. Slicing of graphene shown for two honeycomb cells 
when the contacts are at zigzag edges. Red and green dots 
indicate the two sublattices A and B, the connecting bonds 
of the sites are solid lines. The vertical slices, right side, 
are marked with vertical gray lines, ao represents the lattice 
constant, indicated in blue. 



Here, ^ is the range of the potential and Fj are the lattice 
sites. The concentration of scatterers is given by nimp = 
MmpZ-^tot; where A'tot is the total number of lattice sites. 
We focus on pairwise uncorrelated impurities, (KiV^') = 
{Vn)Sn,n', sacli vanishing on average, (F„) = 0. 

We use as an input parameter the dimensionless corre- 
lation strength Kq, which is defined by the equation for 
the impurity potential correlation function as 



{vimvm = 



2-Ke 



(41) 



In the limit of dilute impurities, summing over all sites 
i,j yields 



K.= 



LW 



Ntot JVtot 



2,,. EE(^(-^)^('^))- (42) 



(Hvo) iVtot i=i j=i 
Inserting Eq. (40) into Eq. (42) and using the relations 

aVslw 



Ntot 
we get^^ 



ao 



and Up = — , 



Til 



Ko = [SVfg, 



(43) 



(44) 



where 



ylf: 



Ni, 



9 \tj TVtotiVimp — 



E 



■A^t. 



E' 



{\^^-R'^\)' 



. (45) 

Thus, the parameter g depends on L, W ^ and ^. It is use- 
ful to introduce the typical impurity strength amplitude 
Vtai which is related to Kq and g by 



v,a^^(y^^^l-{5vY = 



1 




(46) 



VI. CONNECTION BETWEEN THE TWO 
IMPURITY DESCRIPTIONS FROM SECTION 

(III) AND (VC) 

A. Impurity potential 

In this subsection, we establish the link between the 
description of the impurities in graphene in the pseu- 
dospin and isospin representation, as introduced in Sec. 
Ill, and the Gaussian impurities introduced in the previ- 
ous section. ^^'^^ McCann and coworkers assume that the 
different components of the impurity potential Eq. (11) 
are uncorrelated. 



{VrA'?)V,,,,{r')) = VlMr-r')6.J, 



(47) 



where i,j denote the pseudospin and isospin indices as 
introduced in Eq. (11). We can decompose the potential 



due to one impurity at a given site of the lattice, 1^( r^), 
in Fourier components. Defining the vector connecting 
sublattices as to, we find for the Fourier component of 
Vn{^) on sublattice A^'^ 



%,. 



\/3, 



flo 



E^"0 



-ijh~q^ 



and for sublattice B 






V3i 



a-o 



E 



Vni^ — ^)e' 



--i/h~qJ^ 



(48) 



(49) 



where the sum is over all elementary cells. We assume 
V-^^n and V^ „ to be slow functions of the momentum ~q. 
Thus, the quantity Vy.„ is proportional to the first-order 
scattering amplitude for electrons on the same sublattice 
where the impurity resides, while V-^ „ is the one for elec- 
trons on the other sublattice. We will explicitly consider 
two values: the intravalley scattering, q ~ 0, and the in- 
tervalley scattering, g = fcg, where ko connects the two 
different valleys in the reciprocal space and has the am- 
plitude fco = 2/i/3ao. We thus have Vo^n and Vq „ for in- 
travalley scattering in K-K' space and Vko.n and V^ „ for 
intervalley scattering. Intravalley scattering means that 
the electron that is scattered does not leave the original 
cone-shaped valley in k-space, while in the intervalley 
scattering process the electron starts in the K valley and 
ends after scattering in the K' valley or viceversa. 

Equations (48) (49) can now be displayed in the form 
of Eq. (11) using a 4 x 4 matrix notation. For impurities 
Vn{T^) located on sublattice A in the elementary cell r^ 
and small Ij, we can approximate 



q .n 
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(50) 



and, equivalently, for impurities in sublattice B, 
VS = 

q ,n 

















^Cn^" 



-2ik^)r^ 










KnJ 



(51) 



where we used the fact that V^ vanishes due to the 

symmetry of graphene (otherwise we would not have 
any zeros in the secondary diagonal, but terms involv- 
ing V^ ^-^ik'oKy 

As the next step, w convert Eq. (50) and Eq. (51) 
into one single impurity potential matrix that can be 



compared with Eq. (11). This is done by calculating 
the autocorrelation function {{Vif (g) V--^)), with V-^ = 



{{V, ® V^^)) 



(52) 
where the averaging ((•)) is with respect to the positions 
of the impurities and the impurity strength y„, whereas 
the averaging (•) is only with respect to Vn- The nor- 
malization factors are already included in the prefactor 
of ther.h.s. of Eq. (52). 

Comparison of Eq. (52) with the matrix notation of 
the impurity as introduced in Eq. (11) leads, with the 
help of Eq. (46) and the use of the impurity position 
averaged parameter (g) from Eq. (45), to the following 
set of equations 
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(53) 
(54) 
(55) 



(56) 
(57) 
(58) 

(59) 
(60) 



In the limit of short-range impurities, ^ — > 0, we find 



Vo = V^ 



V, 



K = ^4 = o 



(61) 
(62) 



B. Scattering rates 

We now calculate the scattering rates of Eq. (14) and 
Eq. (15) and combine them with the potential Fourier 
components of Eqs. (53), (54), and (55). The density of 
state per spin and per valley j/, which is given by 






(65) 



is inserted into the general expression for the scattering 
rate, Eq. (12), where 



Vq 



(66) 



and Apc = o,q\/3/2 is the area of the primitive unit cell. 
Here we only consider low energies and so we assume that 
Ep ^ Vbg. 

In case of Gaussian scatterers of Eqs. (53), (54), (55), 
and Eqs. (50) and Eq. (51), only the main- and off- 
diagonal matrix elements of Eq. (11) are non zero. This 
helps us simplify Eqs. (14) and (15), so that we find 



1/^3^1 = l/^I^ 







resulting in 
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l/r. 
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(67) 

(68) 
(69) 



Inserting Eqs. (53), (54), and (55) and by using Eqs. 
(12) and (16) leads to 
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(70) 
(71) 
(72) 



Now, using Eq. (65) we find 



while for long-range impurities, ^ — > 00, we find 

Vo = V^ = y/2, (63) 

Vk;-vk = 0- (64) 

In that limit, Eq. (64), only Vb,o is not zero. We can see 
here that, due to symmetry, V^ is always and we are 
left with only one parameter V. This is shown in Fig. 3. 
Also displayed are the potential terms Vo + Vj^', Vq — V^' 
and VjT* which, as we have seen above, are proportional 
to the impurity scattering matrix elements Vq 0, V3 3, and 

Within the simplified picture of this subsection, the 
last result is again consistent with short-range scatterers 
mixing valleys and sublattices while long-range scatterers 
mixing only sublattices. 
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(73) 
(74) 
(75) 
(76) 



Since the dephasing rate l/r^ = in the numerical calcu- 
lations, the low-energy cutoff is provided by the Thouless 
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(a) Effective potentials Vb> Vf^^ and Vq in dependence on ^. (b) Potential terms Vb + Vg', Vb — Vq and Vfe^ in dependence 

on 5. 

Figure 3. Vb, V^ and Vq as set by Eq. (50) and Eq. (51). For short-range impurities, Vo and V^ approach 2V because the 
potential is localized only at sites of one of the sublattices. Wider impurities also cover the neighboring sublattice and Vq start 
to increase and Vb, Vfc,, decrease. Also visible is the faster decrease of Vb — Vq than V^. ^ in units of the lattice constant ao- 



where A is the length or width of the sample, depend- 
ing which one is smaller, ttr = 2to « 2too, if we have 
^/tqo 3> i/Tij, for i,j 7^ 0, and D is the diffusion constant 
given by^^ 
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VqTtb.- 



(78) 



Thus, we get the following ratios: 



r^_VliVo + V^rA'A 
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(79) 
(80) 
(81) 



and Vf*^ depend 



Since the factors (Vb + Vo')^ (Vb -V^f 
on the impurity parameters V and ^, as shown in Fig. 3, 
by tuning one of these parameters and the system size A, 
and gate voltage Vbg, one cam move in the r^/ri - t^/t^ 
diagram of Fig. (1). This phase diagram is equivalent to 
what has been shown in the experimental paper, '^ where 
the temperature has been varied to reach the different 
regions of the diagram. By tuning these parameters, a 
change from weak localization to weak antilocalization 
can be observed. We note that a change of the impurity 
concentration changes t^/ti and t^/tz, while the ratio 
originating from the warping rate t^/t^, is independent 
of the impurity concentration riimp. 

For the magnetic rate 1/tb we find accordingly. 
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TB 



GABeh^v^ 



(Vb + V^)^Apcnin,pVBG ' 



(82) 



VII. NUMERICAL RESULTS 



All results presented in the following are for samples 
with armchair edges. Most calculations have also been 
done with zig zag edges, especially in the case of small 
systems, but they do not show any significant difference 
and are consequently not displayed here. We consider al- 
ways samples with an aspect ratio of 1, L « W, ranging 
from iV = 20 and M = 48; M = 30, N = 72, up to 
M = 80, A^ = 192. The correlation length ^ is given in 
units of flQ, the lattice constant. The correlated disorder 
strength A'o is the dimensionless parameter, defined in 
Eq. (44). For small systems, we use £,/ao = 0.5,..., 3 
and A'o — 0.5, ...,4, while for the larger system, we con- 
sider only ^/ao = 0.5,..., 2, Kq = 0.5,..., 2. Unless ex- 
plicitly mentioned otherwise, the impurity density is set 
to Tijmp = 0.03, meaning that 3% of the atomic sites 
are occupied by impurities. The experiments on weak 
antilocalization"^ indicate indeed that this is a realistic 
concentration of impurities in these graphene samples. 
We also performed the calculations for higher densities 



up to n 



imp 



0.3. 



Impurities are uniformly distributed across the sample. 
For each impurity realization we calculate the electrical 
conductivity as delineated above. We run the numerical 
calculations for Nc = 5000 different realizations in or- 
der to average the results for the conductivity and ther- 
mopower. For the larger systems, Nc = 1000 turns out to 
be sufficient, since self-averaging improves with increas- 
ing system size. The back gate voltage, Vbg is displayed 
in units of the hopping parameter t, which is set to 2.7 eV. 
We consider the range Vbg/^ = 0, ..,0.1, where Vbg = 
corresponds to the Dirac point. The magnetic field is dis- 
played as the magnetic flux $ through the whole sample 
in units of the magnetic flux quantum $0- 
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Figure 4. Average conductance G, as calculated numerically with the recursive Green's function method for dimensions A'^ = 30 
and M — 72. The back gate voltage Vbg is varied as indicated in the figure, (a) Positive magnetoconductance, weak localization 
effect at small correlation length ^ = 1. (b) Negative magnetoconductance, weak antilocalization effect at larger correlation 
length £, — 3. (c) Change between positive and negative magnetoconductance by tuning Vbg- The lines are merely a guide for 
the eye. 
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(b) Size 20 X 48, 5 = 1, K = 2 

Figure 5. Samples of numerical results showing strongly non- 
monotonic magnetoconductance. 



A. Conductivity 

1. Magnetoconductance sign 

Numerical results for the niagnetoconductivity are dis- 
played in Fig. 4. We can distinguish the weak locahza- 
tion from the weak antilocalization by the sign of the 
magnetoconductance: 



ACTsg„(B) = sgn 



fdAcr{B)\ 



Negative magneto conductance Ao'sgn(-S) = — 1 at weak 
magnetic fields B corresponds to weak antilocalization, 
since the magnetic field reduces the quantum correction 
and thereby the conductance. Positive magneto conduc- 
tance A(Tsgn(B) = 1 corresponds to weak localization. 
As seen in Figs. 5, this sign can change at large mag- 
netic fields. This occurs as soon as the magnetic rate 
1/tb exceeds all symmetry breaking rates l/r^ defined 
in Eq. (12). We note that the two-dimensional sam- 
ples considered in the numerical calculations are rather 
small, with length L and width W that are much smaller 
than the magnetic length Is for moderate magnetic fields 
B {Ib = 0.026fim^T/B). Therefore, the sensitivity of 
the conductivity to an external magnetic field B is re- 
duced and the magnetic rate 1/tb becomes suppressed 
for /s > i, VF to^" 



1/ n^^ 

1/tb =cD—^, 



(84) 



V dB 



-y 



(83) 



where c is a geometrical factor of order unity. On the 
other hand, the symmetry breaking rates l/r^ , Eq. (12) 
do not depend on the system dimensions L and W . Only 
the relaxation rate originating from the warping term, 
I/tu,, acquires a small sample size dependence and be- 
comes smaller when Dtw ^ LW . Therefore, as seen in 
Fig. 4, the change of the weak localization corrections oc- 
curs only when the magnetic fluxes through the samples 
corresponds to huge magnetic fields. This is due to the 
small samples considered in the numerical calculations. 

As seen in Fig. 4(a), for small correlation length one 
observes positive magnetoconductance, while at a larger 
correlation lengths negative magnetoconductance occurs 
for the same impurity strength Kq, see Fig. 4(b). As 
we will analyze in detail below, this can be explained by 
the reduction of the intervalley scattering amplitude Vk^ 
with the increase of the correlation length ^, as shown in 
Fig. 3. 

A transition between positive and negative magneto- 
conductance is observed as function of the gate voltage 
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Figure 6. Schematic plots of different types of magnetocon- 
ductance cy{B). At Bmax a-U weak localization corrections are 
suppressed. If the magnetoconductance is monotonous as in 
cases (1), the amplitude is Aa{Bmax)- Otherwise, the mag- 
netoconductance has a maximum at _B = B^x, as in case (2). 



VbGi f-s seen in Fig. 4(c). Similar transitions can be ob- 
served when changing Kq or riimp , which we do not show 
here. 



2. Magnetoconductance amplitude 

In addition to the sign of the magnetoconductance, the 
amplitude of the magnetoconductance Act is important 
and can reveal more information on the nature of the im- 
purities in the sample. As a measure of the amplitude of 
the magnetoconductance, we take the difference between 
the conductance at the first extremum Bg^ and the one 
at zero magnetic field (B=0), 



Aa{Be^)~aiBe:,)-aiB = 0). 



(85) 



Positive A(j(Bex) thus means that there is a weak local- 
ization dip at weak magnetic fields, while negative val- 
ues correspond to the amplitude of the weak antilocal- 
ization peak. We summarize the different types of mag- 
netoconductivity schematically in Fig. 6. The amplitude 
Aa{Bex) is indicated by arrows, both in the case when 
there is no sign change AasgniB) as the magnetic field is 
increased (1), as well as when the sign changes, (2). The 
quantum conductance corrections vanish at large mag- 
netic fields, when Ib is of the order of the elastic mean 
free path l^. This magnetic field we denote as Bmax, as 
defined by Ismax — 'e- 

As seen in Fig. 4, there are still statistical fluctuations 
on smaller magnetic field scales, despite the large number 
of realizations Nc = 5000, which we used in the averaging 
of the conductance. 



Phase 



a. Magnetoconductance sign 

In order to study the crossover between positive and 
negative magnetoconductance as function of the three 
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Figure 7. Numerical magnetoconductance sign phase dia- 
grams calculated with Eq. (83). (a), (b) as function of corre- 
lation length ^ and impurity strength parameter Ko at fixed 
gate voltage, (c) as function of Ko and Vbg for fixed f , and 
(d) as function of ^ and Vbg for fixed Kq. Positive magne- 
toconductance is indicated by the red area, while negative 
magnetoconductance is yellow. The system size is fixed to 
M = 30 and iV = 72. 



impurity parameters ^, A'o, and Vbgi we use Eq. (83) 
at weak magnetic fields to assign the value Acrsg„(i3 — >■ 
0) = — 1 for negative magnetoconductance (NMC) and 
Acrsgn(i? — > 0) = 1 for the positive magnetoconductance 
(PMC). With this information we create three such sign 
phase diagrams, varying two parameters while the third 
parameter stays fixed. For better visualization, the dia- 
grams are colored yellow for negative and red for positive 
magnetoconductance. The crossing line is highlighted by 
the black line. We note that the numerical data is given 
only at the crossings of grid lines. We used an interpo- 
lation method to get the continuous crossing line. As 
seen in the ^ — Kq diagram, PMC occurs for short-range 
impurities with f < an in the whole range of impurity 
strengths Kq . This is expected since short-range impuri- 
ties cause intervalley scattering and thereby suppress the 
pseudospin triplet Cooperons. 

For impurities with larger correlation lengths (^ > oq), 
NMC is observed for strong Kq and large Vbg- Surpris- 
ingly, there is a regime at weak Kq and small Vbg where 
PMC occurs even for large ^. We observed similar be- 
havior for other sizes M and N which we do not display 
here. In the Vbg ~ Kq diagram we see a change from 
PMC to NMC when moving away from the Dirac point. 
At the Dirac point we observe a change at Kq w 1. For 
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0.05 
Vbg 



Figure 8. Vbg — Ko — £, phase diagram for M — 30 and TV = 
72, where the magnetoconductance sign is obtained from the 
numerical results with Eq. (83). The surface with vanishing 
magnetoconductance is displayed. 




V"bg > 0.05 t only NMC is observed for £_ = 2ao. The 
^ — Vbg phase diagram shows that for ^ <C ao PMC oc- 
curs independent of the back gate voltage. At larger ^ a 
change to NMC occurs, as expected. 

In order to get a simpler representation of the results, 
we merge all three two-dimensional phase diagrams into 
one three-dimensional phase diagram, in which each of 
the parameters ^, Kq and T^bg is represented by one axis. 
Fig. (8). 

b. Magnetoconductance Amplitude Using Eq. (85) 
we determine the amplitude of the magnetoconductance 
and plot it in contour plots in the ^ — A'o, the Vbg ~ Kq 
and the ^ — Vbg plane, respectively, in Fig. 9. In these 
figures, the contour lines are tagged with the correspond- 
ing numbers calculated from Eq. (85). The main features 
already observed in the magnetoconductance sign phase 
diagrams can be seen in Fig. 9: PMC occurs for short- 
range impurities with ^ ^ oq irrespective of -ftTo s-nd Vbg- 
For impurities with larger correlation lengths (^ > oq), 
there is NMC for strong Kq and large Vbg, as expected. 
At weak Kq and small Vbg, PMC occurs even for large 
^. Its amplitude remains weak, however, corresponding 
to only a small weak localization dip at weak magnetic 
fields. This is displayed in Fig. 10. 



B. Thermopower 

Next, we calculate the thermopower S, as introduced 
in Sec. IV B, using the Mott formula of Eq. (32) in 
order to see if there are quantum corrections and, as a 
result, magnetothermopower. To this end, we use the 
results for the conductance G and its dependence on gate 
voltage Vbg from the previous section to calculate and 



Figure 9. Magnetoconductance amplitude obtained from the 
numerical results with Eq. (85) in (a), (b) the ^ — A'o for two 
Vbg, (c) the Aq — Vbg and (d) the ^ — Vbg plane, respectively. 
The black line indicates vanishing magnetoconductance. Con- 
tour lines are in units of 2e^ /h. M = 30 and A'^ — 72. 
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Figure 10. Numerical results as function of the parameters 
Vbg, ^0,5 for the magnetoconductance amplitude, Eq. (85), 
for M = 30 and iV = 72. Surfaces of value 0.3 (red), (or- 
ange), —0.3 (yellow) are displayed, in units of 2e^/h. Each 
grid line crossing corresponds to a numerically calculated 
value. 
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Figure 11. Numerical results for (a) dG /(IVbg in units of 
{2e^/h)/eV as function of magnetic flux </> (M ^ 30, N = 72, 
Ko = 2). (b) Magnetothermopower in units of /iV/K^ for 
same system as (a) . The corresponding conductivity is shown 
in Fig. 4(b). 



analyze 



dG 



dVi 



5G 



and the thermopower S. We note that by 



expanding m the quantum corrections of the conductivity 
to first order, we can write the quantum correction to the 
thermopower SS as 






T Ga 



7r2 fc| 6dG/dE 

TM Gel ' 



(86) 



where ^ ..^y = 0.024iJ.V/K^. Since the conductance 

change with magnetic field is only of the order of e'^/h, 
the first term in Eq. (86) is of order ■^. Moreover, 
since Sd vanislres at the Dirac point and increases lin- 
early with the gate voltage, we expect the last term 
to become more important in the vicinity of the Dirac 
point. This is confirmed by the numerical results shown 
in Fig. 11(b) where -;^— is plotted as function of mag- 
netic field. This quantity decays with magnetic field; the 
larger the gate voltage Vbg, the larger is the sensitivity 
to the magnetic field. As shown in Fig. 11(a), the ther- 
mopower shows the same qualitative behavior, namely, a 
reduction of the amplitude of the thermopower when a 
magnetic field is applied. For positive gate voltage the 
thermopower is negative, so that we obtain positive mag- 
netothermopower. 

For the largest samples considered, we also observe an 



initial increase of 



dO 

dVB 



at weak magnetic fields, and a 



peak which becomes enhanced with increasing gate volt- 
age (see Fig. 12), followed by a decay to its classical value 
at larger magnetic fields. Correspondingly, in these sam- 



ples the amplitude of the thermopower (see Fig. 12(c)) 
first is enhanced when the magnetic field is applied and 
then it decreases toward its classical value at larger mag- 
netic fields. 

Since the magnetic field suppresses quantum correc- 
tions, we can obtain the classical value numerically from 
the high-magnetic field limit. Since the thermopower is 
fluctuating even in the range where it should saturate 
towards the classical limit, we needed to average over 
the fluctuations of the data points at different magnetic 
fields. For every set of parameters we thereby calculate 
Smean and usc as a definition for the quantum corrections 
the difference to the thermopower without magnetic field, 



AS* := S{B = 0) - S'„ 



(87) 



In the following, we explore the size and sign of these 
quantum corrections of the thermopower, and how they 
depend on the impurity parameters, such as the strength 
Kq and the correlation length ^, and on the gate voltage. 
In Fig. 13(a) we plot the values of AS* as function of 
these three parameters. As we are particularly interested 
on how these quantum corrections to the thermopower 
are related to the weak localization corrections to the 
conductivity, we also plot the magnetoconductance am- 
plitude (MCA) in Fig. 13(b) in units of 2e^/h. We notice 
that the set of parameters giving the most pronounced 
weak antilocalization seems also to yield the strongest 
thermopower enhancement. We will explain and analyze 
this effect in detail in Sec. VIII D. This behaviour can 
also be seen when looking at the MC amplitude phase 
diagrams as function of Vbg in Fig. 9 since the deriva- 
is the highest when there is a strong increase 



tive 



dG 

dVs^ 



of values in the Vbg direction. This is the case in the 
regime of weak antilocalization, while a smaller change is 
observable in the weak localization regime, see Figs. 9(c) 
and 9(d). 



VIII. ANALYSIS OF NUMERICAL RESULTS 

In order to get a better understanding of these numer- 
ical results for the magnetoconductance and the magne- 
tothermopower, we compare them with the analytical ex- 
pression of Eq. (25) and use the scattering rates l/r^ and 
1/Ti and the effective dephasing rate I/t^ as fitting pa- 
rameters. We also attempt an ab initio calculation where 
we use the input parameters of the numerical calculations 
to calculate directly the different scattering rates l/ry, 
insert these in the analytical expression, and compare the 
result with the numerical ones. 



A. Fitting of the numerical results 

For the fitting procedure, we picked the sample size 
A^ = 30 and M = 72, where we observed a clear tran- 
sition between weak antilocalization and weak localiza- 
tion wlren changing the parameters Kq, ^ and the back 
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Figure 12. Numerical results for (a) conductance (b) dC/dVec in units of (2e^//i)/eV. (c) Magnetothermopower in units of 
liV/K' ((a)-(c): Ko = 0.5, ^ = 2). 




(a) M = 60 and A^ = 144. 



(b) M = 60 and N = 144. 



Figure 13. Numerical results for (a) magnetothermopower 
amplitude, AS/T, Eq. (87), as function of Vbg,^o,^- Sur- 
faces with value —0.05 (red), —0.1 (orange), —0.2 (yellow) 
in units of fj,V/K'^ are displayed, (b) magnetoconductance 
amplitude Act, Eq. (85), as function of the parameters 
Vbg,^o,?- Surfaces of value (red), —0.3 (orange), —0.6 
(yellow), —1 (green) in units of 2e^//i are displayed. (The 
grid corresponds to the numerically calculated values.) 
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gate voltage, see Fig. 9. As mentioned before, since the 
magnetic length exceeds the system size L in the mag- 
netic flux range considered, we have to use the appro- 
priate magnetic scattering rate 1/tb, Eq. (84), which 
is strongly reduced by a geometrical factor compared to 
the two-dimension limit. As discussed above, since the 
numerical results are done at zero temperature, the de- 
phasing rate l/r^ can be substituted by the Thouless 
energy, see Eq. (77). 



1. Fitting with a fixed l/r^ 

In this subsection, the two fitting parameters in Eq. 
(25) are taken to be I/t^ and l/r^, while l/r^ is calcu- 
lated directly from Eq. (77). For the parameter range 
C/ao = 1,1.5,2, Ko = 1,1.5,2 and VbgA = 0.01,0.05, 
we find that Eq. (25) fits well. Within this parameter 
region, we find the transition line between weak antilocal- 



Figure 14. The numerical results for the magnetoconductance 
amplitude AG{<p) = G{(t>)-G{0) for M = 30, iV = 72 together 
with the fitting in terms of l/n, l/r^ in Eq. (25). Red ^ = 1, 
green ^ = 1.5, blue ^ = 2, I/bg = 0.01, Ko = 1. 



ization and weak localization. We concentrate our fitting 
procedure to the validity range of the analytical formula 
at small magnetic fiux, where Is > L, Eq. (84). 

We show the results of the fitting exemplary for some 
parameters in Fig. 14(a), where the magnetoconductance 
is positive. In that case of pronounced weak localization, 
the fitting becomes worse, and the uncertainty in the set 
of l/Ti and 1/tz is larger. 

For each parameter set used in the numerical calcula- 
tions, we get a set of fitted values of the rates t^/u - 
Tip/tz which we plot in the diagrams of Figs. 15(a) and 
15(b) as a symbol. The color of the symbols indicates 
the value of Kq, while the symbol itself depends on the 
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Figure 15. r^/n - t^/t^ diagram obtained from fitting of numerical results. Empty symbols: weak localization. Filled symbols: 
weak antilocalization. Black line: obtained from Aa(B) = in Eq. (30). 



correlation length as shown in the legend. Weak localiza- 
tion corresponds to empty symbols and weak antilocal- 
ization to filled symbols. The range of scattering rates 
is log (t^/tj) ss -4...0 and log (r^/r^) « 1...4. Thus, 
the intervalley scattering rates are typically several mag- 
nitudes smaller than the intravalley scattering rates for 
the range of ^ considered. The black line is obtained from 
A<t{B) = in Eq. (30), indicating the transition between 
weak antilocalization and weak localization. With the ex- 
ception of a few data points, we find good agreement. We 
can improve the agreement by fitting also t^ as shown in 
the next section. 



parameter values is displayed in Fig. 14(b) . As in Sec. 
VIII A 1, we display the resulting I/t^, l/r^, and l/r^ in 
a T^/Ti - t^/tz diagram in Figs. 15(c) and 15(d). 

We find that the results fall in the range log [T^/Ti] k, 
—3 ... and log (t^/tz) « —3 ... 3. The agreement with 
the analytical results, the black solid line in Fig. 15, is 
improved. All numerical results fall on the expected side 
of the weak antilocalization to weak localization transi- 
tion line. 



Ab initio calculation 



2. Fitting 1/r^ 

We repeat the fitting procedure for the same sample 
but now taking l/r^ as a third fitting parameter. We 
focus on the range of parameters £,/ao = 1, 1.5,2, K = 
1,1.5,2 and VBc/t = 0.01,0.05. Close to the transition 
the extra fitting parameter does not change the result 
significantly. The result of the fitting procedure for some 



In this section we will attempt an ab initio calcula- 
tion in the sense that we use the input parameters of 
the numerical calculations to calculate directly the dif- 
ferent scattering rates I/t^j, insert them into the analyt- 
ical formula of Eq. (30) for the weak localization correc- 
tions, and compare the latter to the numerical results. 
For this calculation we choose the following parameter 
ranges. The impurity width f is varied from 0.05 oq to 
1.5 oo in steps of 0.05 ao, the impurity strength V from 
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Figure 16. Magnetoconductance in tiie r^/ri - t^/t^ dia- 
gram. The ratios of scattering rates and dephasing rates 
T(t,/Ti, t^/t* are calculated from the parameters ^ (red dots), 
V (blue squares) and Vbg (green triangles) which are varied 
in the range given in the inset. Arrows indicate the direction 
of increasing parameter. The variation with V and Vbg is 
shown for ^/flo = 0.05 and ^/ao = 1.15, ^ is varied at fixed 
V/eV = 0.5 and Veo/i = 0.04. Dashed line: localization 
transition. White area: PMC. Yellow area: NMC. 



Q.'beV to heV in steps of 0.5 eV (corresponding to Kq 
between approximately 0.05 and 2) and Vbg from 0.01 1 
to 0.1 1 in steps of 0.01 1 (note that the analytical result 
is not valid too close to the Dirac point Vbg = 0). The 
impurity concentration is set to 3% as before, if not ex- 
plicitly mentioned otherwise. 



1. Conductance 

a. Magnetoconductance sign transition 

Using the explicit relations between the scattering 
rates and the parameters ^, V, and Vbg, as derived in 
Sec. VI B, in Fig. 16 we show a t^/ti - t^/t^^ diagram. 

In Fig. 16 we display the magnetoconductance trajec- 
tories in the t^/ti - t^/t^, plane by changing the impurity 
parameters and Vbg- For instance, this is done by chang- 
ing ^ for fixed Vbg = 0-01 1 and V = 5 eV, red dots, and 
by changing V and Vbg blue/green dots, respectively, for 
fixed ^ = 0.05 ao and ^ = 1.5 cq. The impurity concen- 
tration is kept fixed, its dependence inferred from Eqs. 
(79) and (80). 

As expected, increasing the range ^ of the impurities 
induces a change from weak localization to weak antilo- 
calization, see the red arrow. When increasing the im- 
purity strength V one moves toward weak localization. 
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Figure 17. Sign of MC as function of V and ^, corresponding 
to Fig. 16 for two different back gate voltages. Red is PMC, 
yellow is NMC. Compare with the numerical results Fig. 7(a). 
Black line: transition between positive and negative MC. 



see the blue arrow. Increasing the back gate voltage one 
also moves towards weak localization, see the green ar- 
row. The diagram shows that t^Jt^, saturates at fixed 
gate voltage Vbg toward a lower limit since the ratio 
of the warping rate 1/tuj and the dephasing rate, Trf,/Tw 
converges to a finite value, Eq. (81), independently of 
V and ^. At small correlation lengths (^ ^ ao) we find 
that the l/r^ term dominates; a V^ dependence on the 
strength of the impurities and a VgQ dependence on back 
gate voltage emerge, in agreement with Eq. (80). In the 
regime of large correlation lengths (^ ^ ao), we observe 
instead a V^q dependence, which is in agreement with 
the dependence of the warping rate of Eq. (81). For 
comparison, we plot the r^JTi - t^/tz diagram without 
the warping term in Appendix B ( Fig. 26). Next, we 
transform Fig. 16 into the form of the phase diagram of 
Fig. 7(a). The resulting S, — V magnetoconductance sign 
diagram is shown in Fig. 17 for two different back gate 
voltages. 

h. Magnetoconductance Amplitude 

Having investigated the dependence of the sign of the 
magnetoconductivity, let us now study its amplitude us- 
ing Eq. (85) to evaluate Eq. (25) with Ffun given by Eq. 
(26). 

The resulting analytical magnetoconductance ampli- 
tude is shown in Figs. 18 and 19. In Fig. 18(a) the 
amplitude shown in Fig. 18(b) is displayed as function 
of T^JTi and T^JT-t,. Since the warping term dominates 
for large correlation lengths (^ > ao) and t^/tw does 
not depend on ^ and V, regions I and II collapse onto 
a line in Fig. 18(b). We can recognize the localization 
transition that we observed in the sign diagrams of Fig. 
17. In addition, we see that in the range of the local- 
ization transition obtained by the numerical calculation 
in Fig. 8, we observe weak antilocalization with a small 
amplitude when close to the Dirac point for small back 
gate voltages and for weak scatterers (small V, /-fo)- A 
small weak localization amplitude is suppressed already 
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Figure 18. (a): MC amplitude analytically calculated, Eq. 
(85), in units of 2e^/h. Black line: transition between PMC 
and NMC. Positive numbers and red color: weak localization, 
(b) MC amplitude as function of r^/ri and t^Jt^,. 



0.1 




Figure 19. Analytically calculated magnetoconductance am- 
plitude, Eq. (85), as function of ^, V and V^bg- Surfaces with 
value —1 (red), (orange) and 1 (yellow) in units of 2e^ /h 
are displayed. Grid lines are guides to the eye. 



C. Quantum Correction to Conductance at B=0 

The quantum correction to the conductance shown in 
Eq. (23) is plotted in Fig. 20 as function of impurity 
strength V (and, correspondingly, K^) and the correla- 
tion length ^. We note that the transition line between 
a positive and a negative quantum correction at B = 0, 
as indicated by the dashed line, does not coincide with 
the transition from negative to positive magnetoconduc- 
tance which is indicated by the thick black line. There- 
fore, there is a region, denoted as II, where the quantum 
correction to the conductance is positive but the con- 
ductance increases with magnetic field, as one would ex- 
pected for weak localization. Only in region I the positive 
quantum correction coincides with the NMC expected for 
weak antilocalization. In region III, the negative quan- 
tum corrections yields PMC as expected for weak local- 
ization. These regions can be related to the different 
types of magnetoconductance sketched in Fig. 6: apply- 
ing a magnetic field in the region II, we expect a non- 
monotonic magnetoconductane, where the conductance 
first increases, reaching a maximum at i? = B^^, and 
then decays toward the classical conductance, as in case 
(2) of Fig. 6. In the other regions, I and III, the MC is 
monotonic, with PMC in region III and NMC in region 
I. 

In Fig. 20(c) the result for the quantum correction 
shown in Fig. 20(a) is displayed as function of r^p/Ti 
and T^JT^t,. Since the warping term is dominant for large 
correlation lengths ^ > ao, and t^Jt-w does not depend 
on ^ and V ^ the regions I and II collapse onto a line in 
Fig. 20(c). 



D. Quantum Corrections to Thermopower 

In this section we use the analytical theory to study the 
quantum corrections to the thermopower and the result- 
ing magnetic field dependence. In particular, we find out 
in which regime these corrections are large and whether, 
in graphene, they are dominated by the quantum correc- 
tions to the back gate voltage slope of the conductance 
rather than by the weak localization corrections to the 
conductance, as in standard metals. 



1. 



]uantum Corrections at Zero Magnetic Field 



by a weak magnetic field. Thus, in the numerical cal- 
culations it is difficult to identify a region with a weak 
localization amplitude. We typically find by comparison 
that the low-energy cutoff in the numerical calculation, 
is larger than the one obtained from the Thouless energy, 
Eq. (77). 



We first consider the amplitude of the quantum correc- 
tions to the thermopower as defined by Eq. (87), namely, 
as the difference between the value at zero magnetic field 
and the classical thermopower. Thus we need to use the 
weak localization correction to the conductance at zero 
magnetic field, 5a {B = 0), see Eqs. (23) and (25), and 
insert it into the Mott formula. 

Inserting the dependence of the scattering rates on the 
back gate voltage, as given by Eqs. (73), (74). (75), and 
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Figure 20. Analytically calculated quantum correction to the conductance at zero magnetic field, Sa {B — 0) in units of 2e^//i 
for two different back gate voltages, (a) and (b). Black dashed line: transition from positive to negative quantum correction, 
Sa (B — 0) — 0. Black line: Transition from PMC to NMC. (c) Same as (a) but displayed as function of r^/ri and r^/r, . 
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Figure 21. The function i9, Eq. (89) (a) for T4,/t^, = (b) 
for t^/t, — 0. The continuous black line indicates the sign 
change of the function t?. Dashed line: transition between 
PMC and NMC. 



(77), we obtain the slope of the back gate voltage depen- 
dence of the conductance at B = as 



dSa (B = 0) _ e2 1 



9Vbg 



nh Vbg 



^ [Ji 1 ^Z 5 ^W ) ! 



where 



i}iU,t,,t^) = -2 + 



2 ^ 4(1 - 1„,) 
1 + 2U l + t,+t,+ L 



-, (89) 



which depends on the parameter ratios ti = T^JTi and 
tz = Tffi/Tz^ as well as explicitly on the warping rate ratio 

We note that while the quantum corrections to the con- 
ductivity are diverging, in the limit of I/t^ — > 0, which 
corresponds to low temperatures and large system sizes, 
the quantum corrections to the slope of the back gate 
voltage dependence of the conductance, Eq. (88), con- 

2 T 

verge to a finite value of order -^ y — , which depends on 
the scatterings rates as follows: 
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Figure 22. •&. Eq. (89), for two different back gate voltages as 
function of ^ and V . As in Figs. 20(a) and 2Q(b), the black 
dashed line indicates the transition from positive to nega- 
tive quantum correction to the conductance, ba {B = 0) = 0. 
Black line: Transition from PMC to NMC. 



1. When the warping term is negligible, f„, w 0, the 
function i? converges to —2 for large ti, the weak local- 
ization regime, as seen in Fig. 21(a). In the weak antilo- 
calization regime of small intervalley scattering {ti <C 1), 
1? turns positive. 

2. When the warping term dominates, and t^ — > 0, 
the function i? is negative and converges for large inter- 
valley scattering ti 3> 1, the weak localization regime, to 
—6 as seen in Fig. 21(b). Remarkably, i? (and thereby 
dd(j/dVBG) remains for t^ ^- negative for all values of 
ti, even in the regime of weak antilocalization. 

In Fig. 22 d is plotted as a function of ^ and V. For 
small back gate voltage, Vbg = O.Oli, t? is positive in 
the whole region of weak antilocalization (correspond- 
ing to phase I in Figs. 20(a) and 20(b)). This positive 
enhancement of i9 (and therefore positive quantum cor- 
rection to the slope of the back gate voltage dependence 
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of the conductance a,t B = 0, Eq. (88)) matches well 
the numerical results of Sec. VII B). In that regime, the 
intravalley scattering rate l/r^ is expected to dominate 
over the warping rate, and the values of ■& indeed agree 
with those obtained in Fig. 21(a) where t^ = 0. At 
higher back gate voltages. Fig. 22(b), the warping rate 
becomes more important, as tw ^ ^bg increases faster 
than tz with Vbg- Indeed, d remains negative for all val- 
ues of ^ and V , in agreement with Fig. 21(b) where the 
intralayer scattering rate is set to zero, tz = 0. In the 
regime I of weak antilocalization, only a slight increase 
of "d is seen, while it remains negative. 

Now, we are in a position to consider the quantum cor- 
rections to the thermopower (55' as given in Eq. (86). For 
small gate voltages, these corrections are dominated by 
the second term in Eq. (86) , resulting in the weak antilo- 
calization regime (I) in a negative correction 6S/T < of 

order ^ ..^y — 0.024:fiV/K^. In the regime of weak lo- 
calization, the quantum correction to the thermopower 
becomes positive. Since the classical magnetoconduc- 
tance increases with gate voltage, we find that the first 
term in Eq. (86) becomes dominant at large gate volt- 
age, and one recovers SS/Sd ~ —SG/Gd, which is char- 
acteristic of standard metals. To study the competition 
between these two terms in more detail, we use for the 
classical conductivity^^ 
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and obtain Sd by inserting Eq. (90) into the Mott for- 
mula, Eq. (31). With Eq. (24) for Sa and Eq. (88) for 
dSa/dVBG, we can use Eq. (86) to find the dependence 
of 5S/T on the impurity parameter and the back gate 
voltage. We obtain 



(91) 

This result is displayed in Fig. 23. The white area cor- 
responds to the regime where croiff <C Acr(_B = 0), where 
Eq. (86) is no longer valid. 

Comparing the results with Fig. 22 we can see that the 
term due to the quantum corrections to the gate voltage 
slope of the conductance, the •& term in Eq. (91) is dom- 
inant in our parameter range. The quantum corrections 
of the conductance, the 9 term, plays a minor role here. 
However, that term gains importance with higher back 
gate voltage. For small back gate voltage the intravalley 
scattering l/r^ dominance is visible in phase I (WAL), 
where we detect a negative amplitude of the correction 
SS/T, while in phase II and III the correction is positive. 
For higher back gate voltage, the warping term 1/r.u, is 
becoming dominant and we can observe a positive cor- 
rection for the complete phase I and II and even a small 
correction for parts of the phase III. 

To see the connection between the electrical conductiv- 
ity correction and the thermopower correction, we plot da 
versus SS/T in Fig. 24. We focus on the range of weak 





(a) Vbg = O.Olt 



(b) ^BG = O.lt 



Figure 23. Analytically calculated SS/T obtained from Eq. 
(91) in units of /iV/K'^ for two different back gate voltages. 
Black line: transition from PMC to NMC. Dashed line: 5a — 
0. White area: weak localization correction to conductance 
exceeds classical value. 



antilocalization (phases I and II). For small back gate 
voltage. Fig. 24(a), we can see the transition from posi- 
tive to negative thermopower correction dS when increas- 
ing the impurity size ^. We find negative 6S <0 in phase 
I while positive 5S > is seen in phase II. An increase of 
the impurity strength is moving the system toward pos- 
itive SS > 0. This behavior changes at higher back gate 
voltage, see Fig. 24(b), since the warping term becomes 
more important. For short-range impurities, the system 
in phase II does not show a clear relation between 6a 
and SS/T. An increase of the impurity strength tends to 
lower the thermopower correction. With increasing im- 
purity range, when ^ sa ag, we observe an increase of the 
thermopower correction with an increase of the conduc- 
tivity correction near the transition from NMC to PMC. 
For longer ranged impurities, the system is in phase I of 
NMC. In that regime we find good agreement with the 
numerical results and observe a clear relation between an 
increase of 6S/T and 5a. We note that the detailed pa- 
rameter dependence may vary, depending on the value of 
the classical conductance fTDig and the low-energy cutoff 

1/T0. 



2. Magnetothermopower 

Next, we consider how these quantum corrections 
change when applying a magnetic field. We focus first 
on the derivative with respect to the back gate voltage. 
To this end, we can use the expansion at weak magnetic 
fields for Aa{B), Eq. (30), and take its derivative with 
respect to the gate voltage Vbg- For the finite-size sam- 
ples used in the numerical calculations, we can substitute 
l/r^ ~ Ec ~ D / K^ and find that the prefactor in Eq. 
(30) does not depend on Vbg- Thus, we find 
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Figure 24. Analytically calculated SS/T - 5a diagram at zero magnetic field displaying the relation between the quantum 
corrections to thermopower and conductivity corrections for the phases of WAL, I and II, for various impurity parameters ^ 
and V , for two different back gate voltages. Dotted black lines: same impurity strength for different ^. 



where 
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K(ti,tz,tw) in the limits of 1/tu, ~ 0, Fig. 25(a), and 
l/r, = 0, Fig. 25(b). 

For higher magnetic fields we can use the full expres- 
sion and find, 



This is a purely positive magnetothermopower whose 
amplitude increases with system size. We display where 
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and ^E'l is the polygamma function *„ with n = 1, also 
known as trigamma function. Performing the derivative 
shows that the pseudospin-singlet isospin-singlet term 
does not contribute. When the magnetic field reaches 
Smax the magnetothermopower amplitude vanishes as 
does also the magnetoconductivity. We find that the 
thermopower correction 5S/T has a peak for higher mag- 
netic fields when the value for i? = is positive. This is 
in accordance with our numerical finding, see Fig. 12. 



IX. 



ANALYSIS OF EXPERIMENTAL RESULTS 



There have been several reports on magnetocon- 
ductance experiments on single layer graphene. The 
weak localization amplitude is typically of the order of 
2e^ /h.^'^'^"^^^^ Tikhonenko and coworkers showed that it 
is possible to observe weak localization up to a tempera- 
ture of 200 K.^ In their work, they found a clear transition 
between PMC and NMC as the temperature is changed. 
An increase of the temperature decreases the MC ampli- 



I 

tude due to the increasing of the dephasing rate 1/t^{T). 
For example, a change from 5 Kelvin to 14 Kelvin in Ref. 
3 reduces the amplitude at gate voltage Vbg = 40y by 
70%. For this reason, to maximize the amplitude, most 
experiments are done at temperatures below lOK. The 
MC effect is visible in all experiments at magnetic fields 
of up to i? = O.IT. In addition, in all experiments there 
are observable magnetoconductance fluctuations, whose 
amplitudes are however much smaller than the total MC 
amplitude. 

We focus on Ref. 3 at the lowest temperature T = 5K 
since it has the largest amplitude. The scattering rates 
reported by the authors are obtained by fitting the ana- 
lytical theory of Ref.'*. They are found to be log[r0/r*] = 
1.3, 2, 2.4 and log[r0/r,] == -0.2, 1.17, 2 for the three gate 
voltages Vg = 7, 20, 40^. 

Using the dephasing rate l/r^ = 0.1 ps~^ as measured 
in Ref. 3 at T = 5K, identifying Vg = Vbg, and taking 
for the impurity concentration 3%, we find a very good 
agreement between experiments and our results shown 
in Fig. 16. We find the same values for the scattering 
rate ratios t,j,/t^ and r^/ri if we set the impurity pa- 
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Log{r^/r-) 

(a) K{ti,tz,tw) in a t^/tz 
T^/n diagram for I/t™ 



-1 

Log(Ti/T,,) 

(b) K(ti,tz,t-ui) in a t^/t^ ■ 
■ 0. T^/Ti diagram for l/r^ 



Figure 25. i<,{ti,tz,tw) displayed for the two limiting cases 
1/t™ = and 1/r^ = 0. Dashed lines; NMC to PMC tran- 
sition from Fig. 1. Positive numbers indicate positive mag- 
netothermopower (PMT), an increase of 5S/T for small mag- 
netic fields. 



rameters to ^ < 0.5 ag and V w 0.5 eV in the ab initio 
calculations outlined above. Thus, we conclude that the 
analysis of the experimental results with the theory al- 
lows us to make detailed predicitions about the potential 
amplitude and range of the impurities in the graphene 
samples. For the sample of Ref. 3 which was produced 
by mechanical exfoliation of graphite and deposition on 
an oxidized Si wafer, we can conclude that the typical 
range of the impurities in that sample is .J < 0.5ao. If the 
impurity concentration were known from an independent 
measurement, their average strength V could be inferred 
from these magnetoconductance measurements as well. 
Using the measured mobility of He = 12000cm^/(Fs) 
one could in principle estimate one of these parameters. 

Comparing with the analytical results, Fig. 16, we see 
that the experimental results for gate voltages Vq = 7V 
and Vc = 20V are in the regime where the intralayer 
scattering rate l/r^ is dominant while at a gate voltage 
of 40V this rate is superceded by the warping rate l/r^,. 
A higher accuracy of the data analysis can be achieved 
by fitting the experimental data at different gate voltages 
to the same parameters and making use of the analytical 
gate voltage dependence of r^/ri and t^/t^,, as in Eq. 
(81). 

There have been several reports on the measurements 
of thermopower in single layer graphene. ^^*. We found 
above that the magnitude of the quantum corrections 
to thermopower is of the order of l/iV/T^, see Figs. 
13(a) and 24, which is not much smaller than the mea- 
sured values for the classical thermopower of single-layer 
graphene.^"® We note that the dependence of Sa on 
Vbg strongly depends on the gate voltage dependence 
of the dephasing rate l/r^. At low temperatures the de- 
phasing is dominated by the electron-electron scattering, 
yielding'^^ 
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Figure 26. Magnetoconductance in the r^/ri - t^/tz diagram 
for l/tauw ~ 0. The ratios of scattering rates and dephasing 
rates r^/ri, t^/tz are calculated from the parameters ^ (red 
dots), V (blue squares) and Vbg (green triangles) which are 
varied in the range given in the inset. Arrows indicate the 
direction of increasing parameter. The variation with V and 
Vbg is shown for ^/flo = .05, ^ is varied at fixed V/eV — 0.5 
and VBc/i = -04. Dashed line: localization transition. White 
area: PMC. Yellow area: NMC. 



With Ep ~ Vbg and 1/tq ^ Vbg close to the Dirac 
point, one finds that r^ is independent of the gate voltage 
Vbgi ill good agreement with the experimental results."^ 
High magnetic field thermopower measurements at 
room temperatures have been performed. '^^ The tempera- 
ture scale of these experiments ranges from room temper- 
ature down to several Kelvin, which would be the right 
range to observe the magnetothermopower due to quan- 
tum corrections, as studied here, but only when the mea- 
surements would be performed at magnetic fields below 
0.1 Tesla. We are not aware that such measurements 
have been performed to date. 



X. CONCLUSIONS AND DISCUSSION 

We have studied the quantum corrections to the con- 
ductivity and the thermopower in monolayer graphene 
by analyzing numerical results in conjunction with the 
analytical theory. The quantum corrections to the ther- 
mopower result in large magnetothermopower, which we 
demonstrate to be a very sensitive measure of the impuri- 
ties in graphene. While there are experimental measure- 
ments of magnetoconductance of single-layer graphene 
which could be used to determine the average range and 
strength of the impurities, future measurements of the 
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magnetothermopower could provide additional informa- 
tion about the graphene samples. The strong magnetic 
field dependence of the thermopower, the magnetother- 
mopower, is a direct measure of the quantum corrections 
to the thermopower. In contrast to usual metals, its am- 
plitude in graphene is not simply related to the ampli- 
tude of the electrical conductivity, but is also governed 
by the quantum correction to the gate voltage slope of 
the conductivity. 

We demonstrate that the numerical calculations not 
only can be fitted with the analytical theory, but we also 
find even good agreement when the scattering rates en- 
tering the analytical theory are calculated directly from 
the numerical input parameters. This allows a more de- 
tailed understanding, since the dependence on gate volt- 
age and other variable parameters can be studied, reveal- 
ing more accurate information on the impurity parame- 
ters in graphene. 

Besides the sign of the magnetoconductance, showing a 
transition between PMC and NMC, we also analyzed the 
amplitude of the quantum corrections. We found that the 
transition from PMC to NMC does not coincide with the 
transition from positive to negative quantum corrections 
at vanishing magnetic field B — Q. We analyzed the 
relation between the amplitude of the quantum correc- 
tions to the conductance and the magnetothermopower. 
It appeared strongest close to the Dirac point and for 
long-range impurities. 

Here, we studied the thermopower using the Mott for- 
mula, which is only valid for small ratios T/Ep. Still, 
experiments show that it remains valid close to the Dirac 



point. °'^ We intend to study its validity in a future pub- 
lication by using directly the Kubo formula for the ther- 
mopower. For the aspect ratio 1 studied in this paper, 
we did not find a strong sensitivity to the form of the 
boundary conditions, armchair or zigzag. The quantum 
corrections are expected to increase for larger aspect ra- 
tios. Also the warping rate I/tu, is expected to be sensi- 
tive to the aspect ratio, being suppressed when the width 
of the graphene samples is reduced. This effect is similar 
to the reduction of the Dyakonov-Perel spin scattering 
rate in quantum wires. ^^ The other rates l/r^ and r^ are 
expected to depend only weakly on the aspect ratio. We 
note that we assumed in this paper that the dominant 
impurity scattering processes are elastic. The effect of 
magnetic scatterers on the weak localisation corrections 
in graphene will be studied in a separate work. 
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Appendix A: Matrix notation of the Hamiltonian 

For a better understanding of the structure of our for- 
mulation, it useful to write the Hamiltonian in matrix 
notation. We present here Hi , H2 and Himp in this way. 
In particular, the Himp matrix makes it easier to see the 
connection of each part of the impurity potential to the 
isospin-pseudospin description and the matrix notation 
of the Gaussian potential which is shown in section (VI) . 
We have 



Hi ~ Vp(s ® (7 • fc 





= Vp 







U '^x iiK' 
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The second term of the impurity Hamiltonian, as given 
in Eq. (11), can be rewritten in matrix form as 
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Appendix B: Magnetoconductance for XJt^ 







In Fig. 26 we present the results of the ab initio 
calculations for the MC as a function of t^Jti and t^Jt^. 



